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ABSTRACT 

Single crystal ZnO nanowires doped with indium are synthesized via the laser-assisted chemical vapor 
deposition method. The conductivity of the nanowires is measured at low temperatures in magnetic 
fields both perpendicular and parallel to the wire axes. A quantitative fit of our data is obtained, 
consistent with the theory of a quasi-one-dimensional metallic system with quantum corrections due to 
weak localization and electron-electron interactions. The anisotropy of the magneto-conductivity agrees 
with theory. The two quantum corrections are of approximately equal magnitude with respective 
temperature dependences of T~ m and T~ m . The alternative model of quasi-two-dimensional surface 
conductivity is excluded by the absence of oscillations in the magneto-conductivity in parallel magnetic 
fields. 



With the advent of nanotechnology, semiconducting nanowires have attracted extensive interest as the 
building blocks for nanoelectronic devices l ' 4 . ZnO and its doped form have been recognized as a 
promising material system with versatile functionalities based on its remarkable electric, piezoelectric, 
optoelectric and magnetic properties. Undoped ZnO exhibits n-type semiconducting behavior, 
originating from native defects, mainly of Zn interstitials, oxygen vacancies or hydrogen interstitials. It 
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has been found that doping impurities, such as Mg, 5 Ga, 6 In, 7 " 9 P, 10 Sn, 7 and Al, 11 can significantly 
enhance the electrical conductivity, leading to promising applications in nanoelectronic devices. On the 
other hand, nanowires' quasi-one-dimensional (Q1D) structure provides a natural paradigm for studying 
the fundamental physical principles. In particular, the transport conduction theories, such as hopping 4 ' 12 
and weak localization, 13 ' 14 have been reported in both as-grown and doped ZnO nanowires. In addition, 
electron localization in various nanomaterials, such as Si, 15 Sn0 2 , 16 InAs 17 based nanowires, has been 
investigated at low temperatures. Although increasing effort has been devoted to this field, an in-depth 
comprehensive study of the electrical transport is still missing and the underlying mechanisms remain 
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largely unclear. With the aim to understand the microscopic mechanisms, we have performed 
experiments on the temperature and magnetic field dependent conductance of In-doped ZnO (In:ZnO) 
nanowires. Combined with theoretical modeling, we find an excellent quantitative fit of our data in the 
temperature range between 4.2 and 10 K, by treating In:ZnO as a disordered metallic system with 
corrections due to weak localization and electron-electron interaction effects in such Q1D structures. 

The In:ZnO nanowires are synthesized by the laser-assisted chemical vapor deposition (CVD) method. 
Pure zinc powder is placed in the center of the furnace, where the temperature is elevated to 650 °C to 
generate zinc vapor in argon (Ar) under 1 atmosphere pressure. Oxygen (23 ppm) diluted by Ar is kept 
flowing at a rate of 220 seem. Pure indium powder is simultaneously ablated from upstream under 
Nd:YAG (yttrium aluminum garnet) laser pulses. In-doped ZnO nanowires are formed on tin-coated (5 
nm) silicon substrates via the catalytic vapor-liquid-solid process. 19 The energy dispersive X-ray 
(EDX) spectrum (as displayed in the lower left inset of Fig. 1) shows that the indium to zinc atomic ratio 
is around 3%. Selected area electron diffraction patterns (upper right inset of Fig. 1) and high resolution 
transmission electron microscopy (not shown in this paper) reveal the single crystalline structure with 
the (0001) growth direction, indicating that the In atoms are well incorporated into the wurtzite crystal 
lattice. Because of the mismatch of the atomic radii between In and Zn some defects are present. 

The as-synthesized In:ZnO nanowires are suspended in isopropyl alcohol and then deposited onto 
degenerately doped Si substrates that are capped with a 500 nm Si02 layer. An individual nanowire is 
located by scanning electron microscopy (SEM), followed by four-electrode patterning via e-beam 
lithography. Post-thermal annealing at 300 °C for 10 min in vacuum is utilized to further improve the 
contacts. The current-voltage (I-V) and current-gate voltage (I-V g ) curves taken at 4.2 K manifest very 
weak gate dependence. A plot of the resistance versus gate voltage V g is shown in Fig. 1 for two 
measurements. The oscillations persist when a magnetic field is applied but decrease in amplitude when 
the temperature is raised. Similar gate oscillations have been reported in InAs nanowires 17 , attributed to 
the quantization of the momentum of the electrons due to the finite size of the nanowires. However, 
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since the mean free path is found to be much shorter than the radius in our samples, the reason for the 
oscillations is uncertain and needs to be systematically studied. 
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Figure 1 Resistance of a single nanowire versus gate voltage V g at temperature 4.2 K in zero magnetic 
field (two different runs are shown). Lower left inset: the energy dispersive X-ray (EDX) spectrum of 
In:ZnO nanowires. Upper right inset: selected area electron diffraction pattern of an individual In:ZnO 
nanowire. 



The SEM image in the upper left inset of Fig. 2 displays the configuration of a single nanowire (~40 
nm in diameter) in contact with Ti/Au electrodes with 2 \xm inter-electrode spacing for resistance 
measurements. The resistivity of the In-doped nanowire measured by the 4-probe technique at 4.2 K is 
27 liQ m, which is more than two orders of magnitude smaller than that of pure ZnO nanowires. 4 The 
nanowires are in the heavy-doping regime and show degenerate metallic behavior at low temperatures 
with quantum corrections. There are two important quantum corrections: weak localization and 
electron-electron interactions. The key experiments to justify this interpretation are magneto-resistance 
measurements because they allow the derivation of the inelastic scattering length. The magneto- 
resistance is mostly due to a reduction in the weak localization effect, while the two corrections have 
comparable temperature dependence. Using this data, the temperature dependence of the quantum 
corrections is calculated, showing that these corrections account for the temperature dependence of the 
total conductivity. This means that the basic metallic Drude conductivity is constant in this temperature 
range. 
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The magneto-resistance (MR) is measured with the magnetic field ranging from -2.75 to 2.75 T in the 
temperature range from 4.2 to 10 K. Fig. 2 plots the MR data with the magnetic field perpendicular to 
the nanowire axis, and the upper right inset shows the data with the magnetic field parallel to the wire 
axis. As explained in the following analysis section, the negative MR is a consequence of the 
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suppression of quantum interference effects in the weak localization theory. ' It is important to notice 
that there are no oscillations observed in the MR on the scale of 1 T for the parallel field orientation. 
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Figure. 2 Magneto -resistivity (MR) curves obtained with the magnetic field applied perpendicular to a 
nanowire axis at different temperatures ranging from 4.2 to 10 K. Upper left inset: SEM image of a 
single nanowire in contact with four electrodes. Upper right inset: MR plots with the magnetic field 
parallel to the nanowire. 

Figure 3 shows the magnetic field dependence of the conductivity for small perpendicular field from 
to 0.5 T in the temperature range from 4.2 to 10 K. The quadratic slope of the Aa-B 2 curves for the 
magnetic field perpendicular to the wire axis is found to be about twice as large as for that in the parallel 
case (Fig. 3 inset). We will show below that this behavior is in agreement with weak-localization 
theory. To our knowledge, this is the first experimental verification of this predicted asymmetry, since 
the previous experiments focused only on the perpendicular field orientation. 
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Figure 3. Magnetic field dependence of Act = a(B,T)-a(0,T) in the small field range with the 
nanowire perpendicular to the field (for T = 4.2, 6, 8, and 10 K, respectively). The inset shows the 
magneto-conductivity for parallel fields. The quadratic slope for perpendicular field is found to be about 
twice as large as for parallel field. 



Figure 4 plots the conductivity in small magnetic fields ranging from to 0.5 T as a function of 
temperature between 4.2 and 22 K, showing a power law behavior of T~ m . 
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Figure 4. Temperature dependence of conductivity for the field perpendicular to the wire axis with 
values of B = 0, 0.25, and 0.5 T, respectively. The inset shows the data for parallel fields. In the 
temperature range from 4.2 to 22 K the conductivity decreases as a function of temperature with an 
approximate power law of T~ V2 . 



Due to the large surface-to-volume ratio, it has been a question whether the conduction is dominated 
by the surface layer. It was proposed in refs. 13 and 14 that the temperature dependence of the 
conductivity at low temperatures of similar nanowires was proportional to ln(7), and therefore it was 
suggested that these temperature corrections might be due to two-dimensional weak localization in a 
surface layer of their wires. However, our experiments on the magneto-conductivity in parallel fields 
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rule out this interpretation. According to the experiments of Sharvin and Sharvin on hollow metal 
cylinders, if the phase relaxation length is large enough relative to the wire radius, then one would 
expect to see magneto-conductivity oscillations in parallel field with a period of the flux quantum for 
pairs 55* = h 1 2e , where S is the wire cross-sectional area. For our wires of radius 20 nm, the 
corresponding field is B = 1.6 T. As shown in the inset to Figure 2, there is significant magneto- 
conductivity, but no such oscillation effect is observed. Therefore, we instead analyze our results based 
on solid metallic behavior in this temperature range using the theory of quantum corrections to the 
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conductivity due to weak localization and the electron-electron interaction. ' In our case the inelastic 
diffusion length is longer than the wire radius, which sets the system in the quasi-one-dimensional 
regime. 
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The normal-state characteristic parameters for our samples are calculated in the following manner 
using three input parameters. From the data in Figure 4, we assume that the zero-field conductivity in 
this low temperature range is given by the sum of two terms a = cr + Act , where cr = ne 2 r I m* is the 
temperature-independent Drude conductivity and Ac is the temperature-dependent quantum correction. 
By extrapolating the linear behavior of our data to T~ x 12 = we get our first parameter 
cr = 4.25x1 4 SI m Our second parameter, the carrier concentration n, is obtained from the 

dependence of the resistance on the gate voltage shown in Fig. 1 . For a typical doped semiconductor 
one obtains the carrier density n by assuming that the mobility ju = ez I m is constant as the gate voltage 
is changed. However, this assumption is not valid when the carriers form a degenerate Fermi gas. The 

elastic scattering time z in a metal is given by the formula 25 — = — ^( f | u(0 ) f dQ. , where n t is the 

T (27ZJ J 

impurity concentration, u is the Fourier transform of the impurity scattering potential, 6 is the scattering 
angle, and Q is the solid angle. The Fermi momentum p F is proportional to the carrier density to the 

power 1/3. Therefore, the conductivity a is proportional to n 2/3 , and where A^ is 

a R 3 N 

the total number of carriers in the wire. The change in A" is related to the change in the gate voltage by 

dN 

the capacitance C = e—^- of the wire with respect to the gate plate. The usual formula 

C I L = 2tz£ I ln(2/? / r) for the capacitance per length L of an infinitely long thin wire of radius r located 
a distance h above a conducting plane is obtained using the image method. When immersed in a material 
with relative dielectric constant s r , which for SiC>2 has the value of 3.9, the formula for C is multiplied 

by s r . However, in our case, the SiC>2 only fills the region between the wire and the plate and not the 

region on the other side of the wire, which is vacuum. In this case a reasonable approximation would 
be to multiply the above formula for C/L by the average dielectric constant (3.9+l)/2=2.45. Actually, it 
is possible to solve exactly the problem of an infinitely long thin wire above a conducting plane with 
SiC>2 between the wire and the plane by taking the Fourier transform in the direction parallel to the plane 
and solving the resulting one-dimensional equation in the direction perpendicular to the plane. The 
result is that the effective dielectric constant between the wire and the plate is 2.66. Consequently 
C/L = 3SpF/m and C = 16aF for our length L = 2jum . Then from Fig. 1, we use 

dRldV g =45Q/K and i? = 4.3xl0 4 Q and find that Af = 3.0xl0 5 and n = 1.2xl0 26 m" 3 . Of course, the 

charges induced by the gate voltage are confined to a narrow surface layer by Thomas-Fermi screening, 
but since the surface conductivity is in parallel with the bulk conductivity, the actual distribution of 
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charge should not be important. Finally, the third parameter is the effective electron mass, which is 
taken from the literature as the conduction-band mass m* = 0.29m e , 26 

The additional characteristic parameters for our samples are calculated as follows by using these three 
parameters and the free-electron model. The Fermi wave number is k F = {?>7t 2 ri) m = 1.5xl0 9 /n _1 , the 

Fermi momentum is p F = hk F = 1.6xl0~ 25 kg m / s , and the Fermi velocity is 
v F = p F I m* =6.1xl0 5 m/5 . Thus the Fermi energy is E F = 0.3\eV , which is much larger than the 
thermal energy at 10 K, leading to a degenerate Fermi gas, as mentioned above. The elastic scattering 

time is calculated to be r = — — = 3.6 x \0~ l5 s . The mean free path is £ = v f t = 2.2 nm . Notice that the 

ne 

product k F £ = 3A, setting the system in the metallic regime. Finally, the diffusion constant is 

determined to be D = v F £/3 = 4.5 x 1(T 4 m 2 1 s . 

Next we consider the quantum corrections to the conductivity, starting with weak localization. The 
standard formula for the conductivity correction due to weak localization for a wire with a square cross 
section of width a in the absence of a magnetic field, magnetic impurities, and spin-orbit scattering is 

e 2 a ^ 1 

Acr = V j , where S denotes the cross-sectional area of the sample and 

nhS ~ ^(m 2 + n 2 )(7t I of + (L 9 Y 2 



is the inelastic diffusion lengthy = ^Dt p , with diffusion constant D and the inelastic lifetime r p . 
In the one-dimensional regime where a « nL v , the first term in the sum dominates, giving 

e 1 

Acr = L . In the three-dimensional (3D) regime where a»7tL (p , the sum is replaced by an 

nnS 



integral yielding Acr = - 



2n 2 h 



i r 



which depends on a temperature-independent cut-off length i 



on the order of the electron mean free path. In the intermediate range where a ~ jiL 9 , we make a rough 

approximation of replacing the sum by the first term plus the integral, namely by adding the one- 
dimensional and three-dimensional equations. 

The corrections to these equations due to a weak magnetic field have been calculated using first order 
perturbation theory by Altshuler and Aronov. 27 In the above formulas, {L^f 2 must be replaced by 

(L^y 2 +(L ipB y 2 =(L (p0 y 2 +<(2eA/K) 2 >, where A is the vector potential in an appropriate gauge and 

<> denotes the average over the sample. Their calculation is done for a circular sample in a parallel 
field and for a square sample in a perpendicular field in order to easily satisfy the required boundary 
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conditions. The results are that < A 2 >- 



B 2 S 



for the perpendicular field and < A 2 >- 



B 2 S 



for the 



12 Stt 
parallel field. The perpendicular field is predicted to be twice as effective as the parallel field since 

2e 

S7T ~ 25 , which concurs with our measurements. The perturbation theory is valid when — BS < 1 . In 

ti 

fact, for the perpendicular field orientation, there is another correction of order B 2 in second order 
perturbation theory due to the coupling of the vector potential with the component of the fluctuation 
momentum along the wire axis. This correction multiplies the one-dimensional expression for Act by 



1- 



j_ 
30 



2e 

h 



BS 



-1/2 



which clearly indicates the limitation on the validity of perturbation theory 



mentioned above. We find that the three-dimensional corrections to the magneto-conductivity for our 
samples are very small, thus we have used the simple Q1D formula 



Acr(T,B)-Acr(T,0) = 



2 '2e^ 



nhS 2(L B ) 2 nfi 



v " / 
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(L ) 3 to calculate L from the initial slope of the 



quadratic magneto-conductance in perpendicular fields (from Fig. 3). (Note that the conductivity 
corrections due to the magnetic -field dependence of the electron-electron interactions are also 
negligible.) From the magneto-resistance data at temperature 4.2 K, we calculate L -56.3nm. L p0 

decreases smoothly to = 40 J nm at 10 K. If one assumes that L^ is proportional to T p/2 ? the best 
fit to our data gives an exponent p = 0.73 . The corresponding values of the phase relaxation time 

V =0^o) 2 ID are t <po =7.0xl0~ 12 s at 4.2 K and r^ =3.7xl0~ 12 s at 10 K. 

At low temperatures the main contribution to phase relaxation in one dimension comes from electron- 
electron interactions with small energy transfer (compared with k B T), also called the Nyquist noise, 
and the value of the exponent p is predicted to be p = 2/3 , in accord with our measurements. The 



theoretical formula for this Nyquist phase relaxation time is t n = 



h 4 a 2 S 2 



This formula gives 



2e\k B T) 2 D 

t n = 5.6xl0" 12 5 at 4.2 K and t n = 3.1xl0" 12 5 at 10 K. The theoretical values of t n are a factor of 0.8 

smaller than our above experimental numbers for t . 

However, according to ref. 28 the standard theory given above does not apply when the phase 
relaxation is dominated by electron-electron interactions with small energy transfer due to a break-down 
of the quasiparticle description of the Fermi liquid. The revised expression in the one-dimensional case 



isAa(7/,5) = 



2 T 



1.94-^2-^1 



when expanded for small magnetic fields. 30 This formula is 
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the same as the standard theory except for the replacement of L p0 by L N = ^Dt n and the two new 
coefficients: 1.94 and 732. The result of using this corrected formula is that our values of L N are 
reduced by a factor of (32) 1/6 = 1 .78 from the values of L^ stated above, and our experimental values of 
t n are reduced by a factor of (32) 1/3 =3.17 from the values of T p0 given above. Now our experimental 

values are smaller than the theoretical results by a factor of 0.4. It should be noted that in ref. 29, in 
spite of the good agreement of the temperature exponent of t n over several decades, there are similar 
difficulties in fitting the prefactor. 

The revised formula gives almost exactly the same zero-field temperature dependence of Act because 
the two correction factors almost exactly compensate each other, i.e. 1.94/1.78 = 1.09. The magnitude 
of the one-dimensional conductivity correction due to weak localization remains the same, and the 
temperature dependence remains proportional to 3 . The one-dimensional weak localization theory 

2 

gives a conductivity difference between 10 and 4.2K of \L {A.2K) - L (1 0K)~\ = 962 (fim)"' . The 

three-dimensional weak localization theory gives a conductivity difference of 



27T 2 k 



1 



1 



= 84 (Q.m) 1 . Thus, the Q1D effect is much more significant than the 3D 



L^IOK) L 9 {A.2K) 
correction. 

Furthermore, there are still the electron-electron interaction effects to consider. The contributions to 
the conductivity due to the electron-electron interaction (also called the Coulomb anomaly) are similar 
to those due to weak localization but with the length L f replaced by the thermal diffusion length 



hD 



L T = , — , and the coefficients are slightly modified. The exponent of the temperature dependence L T 
V kT 



is -1/2, the same as we have used in Fig. 4. In the Q1D regime, Aa ee = ■ 



0.39 



3F^ 



TThS 



L T , while 



in the 3D regime Act = - 
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From our material parameters we calculate 



the correction factor to be F — 0.6 . Using our estimated value for the diffusion constant 
D = 4.5xl0" 4 m 2 Is , we calculate L T = 28.6«m at 4.2 K and L T = l$.5nm at 10 K. Then the calculated 

conductivity differences between 10 and 4.2 K due to the electron-electron interaction are 751(Qm) _1 
(Q1D term) and 47(Qm)~'(3D contribution). 
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Altogether considering the 3D and Q1D weak localization and electron-electron interaction 
corrections, the sum of the four calculated contributions to Act between 10 and 4.2 K is 1844(Qm)~'. 
The magnitudes of the weak localization and the electron-electron corrections are approximately equal. 
Experimentally our conductivity at 4.2 K is 3.723 xl0 4 (Qm) _1 and at 10K is 3.908 xl0 4 (Qm) _1 , so the 

difference is 1850(Qm)~'. Therefore, we believe that the combined theory is in excellent agreement 
with our data. 

In conclusion, we have presented low temperature magneto-conductance measurements on ZnO 
nanowires highly doped with In. The temperature and magnetic field dependence of the conductivity of 
our wires is in good agreement with the theory of quasi-one-dimensional weak localization and electron- 
electron interactions in a disordered metallic system. In particular, the quadratic magneto-conductance 
slope for the magnetic field perpendicular to the wire axis is found to be about twice as large as for that 
in parallel case. To our knowledge, this is the first experimental verification of this predicted 
asymmetry. In addition, the absence of oscillations in the magneto-conductance curves for magnetic 
fields parallel to the nanowires excludes the previous suggestions and indicates that the conductance of 
these nanowires is a bulk, rather than a surface effect. 
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